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for Mobile Autonomous Robots 
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. Abstract 

If a smooth, closed, and embedded curve is deformed along its normal vector field at a rate 
, proportional to its curvature, it shrinks to a circular point. This curve evolution is called Euclidean curve 

shortening and the result is known as the Gage-Hamilton-Grayson Theorem. Motivated by the rendezvous 
problem for mobile autonomous robots, we address the problem of creating a polygon shortening flow. 
A linear scheme is proposed that exhibits several analogues to Euclidean curve shortening: The polygon 
shrinks to an elliptical point, convex polygons remain convex, and the perimeter of the polygon is 
\Q ' monotonically decreasing. 

I. Introduction 

This paper studies the rendezvous problem for mobile autonomous robots, in which the goal is to 
develop a local control strategy that will drive each robots's state (usually its position) to a common 
value. Research on this problem has been performed in discrete and continuous time. The discrete time 
research can be split further into synchronous systems [l]-[5] (i.e., each robot moves only at global clock 
ticks), and asynchronous systems [6], [7] (i.e., no global clock is present). In the synchronous case there 
have been several papers on circumcenter algorithms [l]-[4], in which each robot moves towards the 
center of the smallest circle containing itself and every robot it sees. In both the continuous and discrete 
time cases, the research has assumed fixed communication topologies — the sensors are omnidirectional 
and have a range larger than their environment, allowing each robot to see all others — and time- varying or 
state-dependent communication topologies — the sensors have limited range; the sensors are directional; 
or, communication links may be dropped or added. In continuous time, a fair amount of research has 
been based on a simple strategy called cyclic pursuit [8]— [1 1] . In this strategy the agents are numbered 
from 1 to n, and each agent pursues the next with the nth agent pursuing the 1st. 

In this paper we look at the rendezvous problem from a different perspective. We are concerned with 
the shape of the formation of robots as they converge to their meeting point. We would like the formation 
to become more "organized," in some sense, as time evolves. We use a simple model, numbering the 
robots from 1 to n and considering a fixed communication topology in continuous time. We then view the 
robot's positions as the vertices of a polygon, and, motivated by the Gage-Hamilton-Grayson Theorem 
described below, we seek to create an analogous polygon shortening flow. 

To introduce the Gage-Hamilton-Grayson Theorem, consider a smooth, closed curve x(p, t) evolving 
in time: p G [0, 1] parameterizes the curve; t > is time; and x(p, t) G M 2 . We can evolve this curve 
along its inner normal vector field N(p, t) at a rate proportional to its curvature k(p,t) (curvature is the 
inverse of the radius of the largest tangent circle to the curve at x(p, t), on the concave side): 

— (p,t) = k(p,t)N(p,t). (1) 

This curve evolution is known as the Euclidean curve shortening flow [12], and is depicted in Fig. [2 
Let L(t) and A(t) denote respectively the length and enclosed area of the curve at time t. Gage [13]- 
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[15], Hamilton [15], and Grayson [16], [17] showed that a smooth, closed and embedded curve evolving 
according to (0 remains embedded and shrinks to a circular point. The term "circular point" means that 
the curve collapses to a point and, if we zoom in on the curve as it is collapsing, the curve is becoming 
circular. Throughout the evolution, A(t) = —2ir and L(t) is monotonically decreasing. In [17] it is also 
stated that under Q, "the curve is shrinking as fast as it can using only local information." This notion 
will be clarified later. 

There has been prior work in creating polygon shortening flows. Motivated by the curve shortening 
theory, and applications in computer vision, Bruckstein et al. [8] study the evolution of planar polygons 
in discrete time. A scheme is proposed that shrinks polygons to elliptical points (the vertices collapse 
to a point, and if we zoom in on the collapsing polygon, the vertices are converging to an ellipse). In 
addition, [8] discusses a polygon shortening scheme based on the Menger-Melnikov curvature [18]. In 
[19] this scheme is studied and it is shown that most quadrilaterals shrink to circular points. In [20] a 
flow is formulated such that the area enclosed by the polygon shrinks at a rate of 2ir and the perimeter 
of the polygon is monotonically decreasing. 

In this paper we study a planar polygon, with vertices z\, . . . , z n in the complex plane C, as it evolves 
according to 

Zi = -z{zi + i-Zi) + -(zi-i-Zi), i = l,...,n, (2) 

where the indices are evaluated modulo n. Thus, vertex i pursues the centroid (center of mass) of its two 
neighboring (according to numbering) vertices. A discrete-time version of © is studied in [8], and it is 
shown that the polygon shrinks to an elliptical point. The contributions of this paper are as follows. We 
introduce the curve shortening theory and its relation to the rendezvous problem. We also demonstrate 
the importance of studying the shape of the formation of robots as they rendezvous. We then show the 
following under ©: 1) convex polygons remain convex, 2) if vertices are arranged in a star formation 
about their centroid, they remain in a star formation for all time (in particular, the robots will not collide), 
3) the perimeter of the polygon monotonically decreases to zero. Finally, we derive the optimal direction 
for shortening the perimeter of a polygon. 

II. Polygon shortening 

We consider n robots in the plane to be the vertices of an n-sided polygon. In this section we formally 
define a polygon and introduce two polygon shortening schemes. 

A. Definition of an n-gon 

Following [21] we introduce the definitions of a polygon and a simple polygon in M 2 (or equivalently 
C). An n-gon (n-sided polygon) is a (possibly intersecting) circuit of n line segments z\Zi, Z2Z3, ■ ■ ■ , z n z\, 
joining consecutive pairs of n distinct points z%, Z2, ■ ■ ■ , z n . The segments are called sides and the points 
are called vertices. A simple n-gon is one that is non-self-intersecting. We denote the counterclockwise 
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Fig. 2. The circumcenter for three points on the curve x(p). 



internal angle between consecutive sides ziZi+i and Zi-\Z\ of an n-gon as (as always, indices are 
modulo n). For a simple n-gon these angles satisfy Ya=i Pi = ( n ~ ^) 7T - An n-gon is convex {strictly 
convex) if it is simple and its internal angles all satisfy < < it (0 < /?, < it). 

B. Shortening by Menger-Melnikov curvature 

We now briefly describe the polygon shortening scheme studied in [8], [19], and our reasons for not 
following this approach. Let x(p), p € [0,1], be a smooth curve. Consider a set of parameter values 
pi < p2 < ■ ■ ■ < p n and the corresponding discrete points x(j>j). By connecting these points we create 
an n-gon. As n — > oo and if the parameter values {pi} become dense in [0, 1], the n-gon converges to 
the smooth curve x(p). The idea is to create a polygon shortening scheme so that as n — > oo, the scheme 
tends to {TJ. 

If three consecutive points x(j>j_i), x(j>j), x(j>j + i) are not collinear, there exists a unique circle (the 
circumcircle) that passes through them. Denote the radius of the circle by R(pi) and the center of this 
circle by C(pi), as shown in Fig. [2] The quantity l/R(pi) is called the Menger-Melnikov curvature and 
has the property that 

1 

lim "757 — V = 

In addition, as the points x(»;_i) and x(pj+i) approach x(j>j), the quantity (C(pi) — x(j>j)) / R{pi) 
approaches N(pj) if k(pi) > and — N(p,) if k(pi) < 0. Therefore, we have 

C(pj) - jc jpi) 

The Menger-Melnikov flow is then given by 

■ / s C{Pi) - x(pj) . 

This flow was studied in [8], [19]. However, due to the complexity of the system the results are 
quite limited [8]. In [19] it is shown that a simple n-gon collapses to a point in finite time, and for 
n = 4 most quadrilaterals tend to regular polygons. However when n is small, this flow may yield 
a poor approximation of the inner normal vector, as shown in Fig. [3] In fact, for a convex n-gon, the 
approximation to the normal vector may not even point into the interior of the n-gon. Also, as the polygon 
collapses, the velocities of the vertices approach infinity, which is not ideal for our application. In light 
of these remarks, we propose the scheme presented next. 

C. Linear scheme 

The linear polygon shortening scheme is given by ©. Defining the aggregate state z = (z\, . . . ,z n ), 
we get the simple form i = Az. The matrix A is circulant (see Davis [22]) and is given by A = 
circ(-l,i,0, ...,0, i), with 



lim 571 \2 = ^(Pi) N (Pi)- 



circ(ai,a 2 , . . . ,a n ) :- 



ai a-2 ■ ■ ■ a n 
a n ai • • • a ra _i 

S2 «3 ••' ai 
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Fig. 3. The normal vector, and the Menger-Melnikov approximation, when the number of points n is small. The approximation 
to N(pi) is very rough. 



The matrix A can be written in terms of the polynomial 

QA{s) = ^s n ~ x + -a - 1, 
and the matrix P = circ(0, 1, 0, . . . , 0), as A = qji(P). By the Spectral Mapping Theorem we obtain 

eigs{A) = {5a(1), <?a(^ 2 ), • • • , gAK" 1 )}, 

where u = e 27rJ//n . Therefore, denoting Aj := g^(cj l_1 ), we have eigs(A) = {Aj : « = l,...,n}. 
Evaluating Aj we get 

A = I e 27ri(n-l)(i-l)/n , 1 2-Kj(i-l)/n _ 1 
* 2 2 
= cos(27r(i — l)/n) — 1, 

where i = 1, . . . , n. From this, one can easily verify the following properties: 1) the eigenvalues of A are 
real, with one eigenvalue at zero, and all others on the negative real line, 2) the centroid z := Y2i=i z i/ n ^ 
stationary throughout the evolution, and 3) the robots asymptotically converge to this stationary centroid. 

The following theorem characterizes the geometrical shape of the points Zi{t) as they converge to their 
centroid and is proved for discrete time in [8], and for general circulant pursuit in [23]. 

Theorem 1: Consider n points, z\{t),... ,z n {t) evolving according to ©. As t — > oo these points 
converge to an ellipse. That is, z\(t), . . . , z n (t) collapse to an elliptical point. 

III. Star formations stay star formations 

We will now see that if a group of robots have the initial arrangement shown in Fig. |5J called a star 
formation, then under © it stays in a star formation and, in particular, no collisions occur. We require 
some preliminary tools. For z G C, we let 9ft{z}, Q{z} and z, denote the real part, imaginary part, and 
complex conjugate of z respectively. 

Lemma 2 (Lin et al. [10]): Let Z\, Z2, and z% be three points in the complex plane, as shown in Fig. 
0] Let n := \z\ - z 2 \, r 2 := \z 3 - z 2 \ and 

F = Z{( Zl -z 2 )(z 3 - z 2 )}. 

Then (i) < a < it, n > 0, and r 2 > if and only if F > 0; (ii) it < a < 2tt, n > 0, and r 2 > if 
and only if F < 0; (iii) the points are collinear if and only if F = 0. 
Proof: Introduce the polar form 

Z\ - z 2 = rie> dl , z 3 - z 2 = r 2 e> 62 
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Fig. 4. The setup for the definition of the function F. 




Fig. 5. A counterclockwise star formation. 



where 0\, 02 are the angles of the line segments in the global coordinate system. Then 

F = Q{(zi — Z2)(zz — Z2)} = ^s{rie~^ dl r2e' e2 } = rir2sm(a). 

Thus, < a < ir, r\ > 0, and r 2 > iff F > 0; and ir < a < 2ir, ri > 0, and r 2 > iff F < 0. Also, 
the points are collinear iff F = 0. ■ 
Now consider our system of n robots, whose positions, not all collinear, are denoted by z\, . . . , z n . Let z 
be the centroid and r.j be the distance from the centroid to z: t . Let ai denote the counterclockwise angle 
from zzi to zzi + \ for i = 1, . . . , n, modulo n. 

Definition 3 (Lin et al. [10]): The n points are arranged in a counterclockwise star formation if > 
and Oj > 0, for all i = 1, . . . ,n, and Y17=i a i = ^ 7r - They are said to be arranged in a clockwise star 
formation if r« > and a-i < 0, for all i = 1, . . . , n, and Ya=1 a i = ~^ 7T - 

This formation is shown in Fig. [5] In what follows we will consider only counterclockwise star 
formations, since the treatment for clockwise star formations is analogous. Also, the case n = 2 is 
trivial, so it is omitted. 

We are now ready to state the main theorem of this section. 

Theorem 4: Suppose that n distinct points, with n > 2, are initially arranged in a counterclockwise star 
formation. If these points evolve according to © they will remain in a counterclockwise star formation 
for all time. 

The proof uses these two results: 

Lemma 5 (Lin et al. [10]): Suppose that n distinct points, z\, ■ ■ ■ , z n , with n > 2 are in a counter- 
clockwise star formation. Then «j < n, Vi 

Lemma 6 (Lin et al. [10]): If n points, Z\, . . . , z n evolving according to © are collinear at some time 
t\, then they are collinear for all t < t\ and t > t\. 

Proof of Theorem |?J We begin by considering the function 

Fi(t) = %{(zi(t) - z)(z i+ x(t) - z)} = r { r i+1 sin(ai). 

By the definition of a counterclockwise star formation we have rj(0) > and < «i(0) < n, Hence 
by Lemma |2j Fi(0) > 0, Vi. We want to show that Fi(t) > 0, Vi and Vt, which by Lemma |2] shows that 
the vertices are in a counterclockwise star formation for all time. 
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Suppose by way of contradiction that t\ is the first time that some Fi becomes zero. We can select 
i = m such that F m (t\) = and F m +i(tx) > 0, for if all the Fj's are zero at tx, then the points are 
collinear, which by Lemma |6] is a contradiction. Hence, we have F{(t) > for all t € [0, t\) and all i, 
F m (t 1 ) = 0, andF m+1 (t!) >0. _ _ _ 

Taking the time derivative of F m , and noting that z = 0, we have F m = ^s{z m {z m+ i — z) + 

{Zm — z)z m+ i}. 

By adding and subtracting z in each term in © we can write © as 

Zi = ~ji Zi+1 ~ ^ + \^ Zi ~ X ~ ^ + ~ 

Using this expression for z m and i m +i and simplifying, we obtain F m = —2F m + G m , where 

G m = ^^{(^m-l ~~ z)( z m+l ~ z ) + ( z m ~ z){z m +2 ~ z)} 

= 2 ( r m-i r m+i sin(a m _i + a m ) + r m r m+2 sin(a m + a m +i)) • (3) 

Now, if F m (ti) = 0, by Lemma one of the following four conditions must hold: (i) a m {t\) = it and 
r m (h),r m +i(h) > 0; (ii) a m {ti) = and r m (ti),r m+1 (ti) > 0; (iii) r m (t x ) = 0; (iv) r m+ i(ti) = 0. 

Condition (iv) cannot hold since F m+ i(ti) > 0. Condition (i) cannot hold, for if it did, all points 
would lie on, or to one side of, the line formed by z m+ \ and z m , a contradiction by either Lemma [5] or 
|6] Assume that condition (ii) holds. Then a m {t\) = and from © we obtain 

G m (ti) = ^ (r m _ir m+ i sin(a m _i) + r m r m+2 sin(a m+ i)) 
_ f m+1 F m _i(ti) H — F m+ i{ti 

^ \ I'm 

Since r m (t{), r TO +i(ti) > 0, F m+ i(ti) > 0, and i 7 ' m -i(ti) > 0, it follows that G m (t{) > 0. By continuity 
of G rn there exists < to < *i such that G m (t) > for all t G [io>*i]- Also, by assumption, F m (t) > 
for t € [0, ti). Therefore F m (t) = —2F m + G m > —2F m for all £ G [to,ti). Integrating this and using 
the continuity of F m , we obtain F m (ti) > e ~ 2 ( tl_t °) F m (t Q ) > 0, a contradiction. 

Finally, suppose condition (iii) holds and r m (t\) = 0. Then z m (t\) is positioned at the centroid, z. 
Assume without loss of generality that 2 = 0. Notice that if Zj(ii) = 0, the angle 6i(t\) is not defined. 
We now establish that if Zi(t\) = and Zi{t\) ^ 0, then lim t | tl 6i(t) is well defined. Expanding z% about 
t\ we have Zi(ti) = Zi(h - h) + hzi(ti) + 0(h 2 ), where 0(h 2 )/h -> as /i -> 0. If ^(ii) = then 
Zi{t\ — h) = —hzi(ti) + 0{h 2 ). Hence, lim^o — = — ij(ti). Therefore the limiting motion 
of Zi(t) as 1 1 ti is along the ray defined by — ij(ti), as shown in Fig. [6] Because of this, we can define 

i(h) if nih) > o, 

- < arcta „ (|tiiM) if n(tl) _ o. <4) 



With this definition we can talk about 6i(ti), and ai{t\), when rj(ii) = 0. 

Suppose that by a rotation of the coordinate system, if necessary, the vector z m+ i{t\) + z m -i(t\) lies 
on the negative real axis. Then we can write 

Z m+ l{tl) + Z m -l{tl) 

— — — - = — r, where r > 0. (5) 

We have r > for if r = then z m _i(ii), z m (ti), 2 m +i(ii) all lie on a line through the centroid, and 
all other points must lie either on or to one side of this line, implying that is not the centroid, or all the 
points are collinear, both contradictions. Since z m (ti) = 0, from (EJ and (J5jl we have z m (ti) = —r, as 
shown in Fig.0 If n = 3 then z m (t\) = and the centroid of z m+ i{t\) and z m -x{t\) is at — r, implying 
that is not the centroid of the three points — a contradiction. 

Therefore we need only consider n > 3. Since z m (t\) = — r, from (@J we obtain 

M*i) = 0- (6) 
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Fig. 6. The limiting 9i(i) as t | ti 
when Zj(ti) = 0. 



Fig. 7. The position of the points z m -i, Zn 
and z m +i at £ = ti, 
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Fig. 8. The required geometry such that m -l(ti) £ [— tt, 0], 6 m +i(ti) £ [0, 7r], and z m +i(ti) + z m -i(ti) — — 2r. All points 
lie either on or to one side of the dotted line. 



To obtain a contradiction for n > 3 we will show that © and © cannot both be satisfied. To do this 
we consider two cases, r m _i(ii) = and r m _i(ti) > 0. Since the points are in a star formation until 
t\, we know that Vi, Oi(t) G (0,tt) for t G [0,t\). Hence, if 6i{t\) and are defined via ©, then 

by continuity, Oj(ti) G [0, tt]. 

If r m _i(ii) = then from $5^ we have <2 m +i(ti) = — 2r. Therefore 9 m+ i{ti) = tt and from ©, 
9rn(ti) = 0. However this implies that all other #i(ii)'s that are defined must lie in [— tt, 0]. Hence 
^{^(ii)} < Vi, which implies that all points are collinear, or that is not the centroid, both 
contradictions. 

If r m _i(ti) > then from and since a m (t\), a m _i(ii) G [0,tt], we have that 9 m+ i(ti) G [0,7r] 
and f? m _i(ii) G [— tt, 0]. So S{z m+ i(ti)} > and 9{,z m _i < 0. Because of this, as can be 
verified in Fig. [8] for (J5jl to be satisfied either z m -i(ii) and z m +i(ii) are both real, in which case 
(ti) = tt, or neither is real and 9 m +i(h) — 8 m -i{ti) > But this implies that all 
points lie on, or to one side of, the line formed by z m -\{t\). Thus all points are collinear, or is not 
the centroid, both contradictions. □ 

Fig. [9] shows the evolution of a polygon that is in a star formation about its centroid. Notice that the 
polygon remains in a star formation, becomes convex, and collapses to an elliptic point. 

IV. Convex stays convex 

We now turn to the case where the formation is initially a convex n-gon. We will show that a convex 
n-gon evolving according to © remains convex. To do this we require a function similar to that in 
Lemma |2 but which measures the counterclockwise internal angle between two sides of an n-gon. 

Lemma 7: Consider a simple n-gon lying in the complex plane, whose vertices Zi are numbered 
counterclockwise around the n-gon. Let /3 2 denote the counterclockwise angle from the side Z2Z3 to the 
side z\Z2 as shown in Fig. and define p\ = \z\ — z 2 \, p 2 = \%3 — 2b\ and 



H = %{(z 1 -z 2 )(z 3 -z 2 )}. 



Fig. 9. The evolution of a polygon whose vertices are in a star formation about their centroid *. The dashed lines show the 
trajectories of each vertex. 




Then (i) < (3 2 < tt, pi > 0, and p 2 > if and only if H > 0. (ii) it < (3 2 < 2ir, p 1 > 0, and p 2 > if 
and only if H < 0. (iii) the points are collinear if and only if H = 0. 
Proof: We introduce the polar form: 

z\ - z 2 = pie' 11 , z 3 - z 2 = p 2 e jl2 , 

where 71, 72 are the angles shown in Fig. EH Then 

H = 3{(zi - z 2 )(z 3 - z 2 )} = 9{/>ie J71 p 2 e~ n2 } = pip 2 sin(/3 2 ) 

Thus, < < pi > 0, and p 2 > iff H > 0; and vr < f3 2 < 2vr, p x > 0, and p 2 > iff H < 0. 
Also, the points are collinear iff H = 0. ■ 

Lemma 8: If an n-gon is strictly convex, with its vertices Zi, i = 1, . . . , n, numbered counterclockwise 
around the n-gon, then these vertices are in a counterclockwise star formation about their centroid. 

Proof: The centroid, z, of the n vertices must lie in the interior of the strictly convex n-gon for if 
it lies on the boundary or in the exterior, we could draw a separating line through the centroid for which 
all vertices lie on, or to one side, contradicting the position of the centroid. With this observation, and 
the fact that the n-gon is convex and numbered counterclockwise, the result is straightforward. ■ 

With these two lemmas, and Theorem 0J we can prove the main result of this section. The result is 
analogous to convex curves remaining convex under (Q, which is shown in [13]. 

Theorem 9: Consider a strictly convex n-gon at time t = 0, whose vertices zi, i = l,...,n, are 
numbered counterclockwise. If these vertices evolve according to ©, the n-gon will remain strictly 
convex for all time. 
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Proof: We begin by considering the function 



Hi(t) = 9f{(zi_i(t) - Zi(t))(z i+1 (t) - Zi(t))} = pi-ipi sin(A) 

By the definition of a strictly convex n-gon we have that pi(0) > 0, and < /3i(0) < n, Vi. Hence by 
Lemma Hi(0) > 0, Vi. We want to show that Hi(t) > for all t, which by Lemma shows that the 
n-gon remains strictly convex for all time. 

Suppose by way of contradiction that t\ is the first time that an Hi becomes zero. We can select i = m 
such that H m (ti) = and H m+ i(ti) > 0, for if all the H^s are zero at t\, then the points are collinear, 
which by Lemma |6] is a contradiction since the points started in a convex n-gon formation. Hence, we 
have Hi(t) > 0, V t e [0,h) and i = 1, . . . ,n, H m (ti) = 0, and H m+1 (ti) > 0. 

Taking the derivative of H m along the trajectories of (O, we have 



Substituting in © for z m _i, i m , z m+ \, and simplifying we obtain ij m = —2H m + G m , where 



1. 



G m — — 9{(^m-2 — z m-l)(z m +l — z m) + ( z m-l — 2 m)(^m+2 — z m+l)}- (7) 

Now, if H m (ti) = 0, by Lemma0 one of the following four conditions must be satisfied, (i) f3 m (ti) = it 
and p m -i(ti),p m (h) > 0. (ii) /3 m (ii) = and p m -i{h), Pm(h) > 0. (hi) p m (h) = 0. (iv) p m _i(*i) = 
0. 

Condition (iii), in which p m {ti) = 0, cannot be satisfied since H m+ i{t\) > 0. Also, since the n-gon is 
initially convex, by Lemma|8]it is in a counterclockwise star formation. By Theorem|4]the vertices remain 
in a star formation for all time and thus remain distinct. Therefore, condition (iv) in which p m -i(ti) = 0, 
cannot be satisfied. 

Assume condition (i) is satisfied. Then (3 m (ti) = vr, /3 m _i(ti), p m {t\) > and H m {t\) = 0, H m+ i(t\) > 
0. Since (3 m {ti) = tt> we have that 

z m +i(ti) - z m (h) Zm—1 (h) - z m (ti) 

Pm Pm—l 

Combining this with the expression for G m we have 

G m (tl) — (z m — 2 Z m —\){z m —\ Zm) (Zfn+i Z m )( K Z rrl J r 2 ^m+l)} 

^ Pm—l Pm 

' ( / '-- J ff m _ 1 (t 1 ) + ^i J ff m+ i(t 1 )V (8) 



^ \Pm—l Pm / 

Since p m ~ x (^i), Pm(h) > 0, H m+1 (t x ) > 0, and iJ m _ x (ti) > 0, it follows that G m (h) > 0. By 
continuity of G m there exists < to < h such that G m (t) > Vt € [io,ii]- Also, by assumption, 
H m (t) > for t G [0,ii). Therefore H m (t) = -2H m + G m > -2H m , for all t <E \to,tx). Integrating 
this and using the continuity of H m , we obtain H m (ti) > 0, which is a contradiction. 

Finally, assume condition (ii) is satisfied. Then (5 m (ti) = and p rn ^i,p m > 0. The angle (3 m is the 
interior angle between the edges z m -\z m and z m z m+ i. For all t £ [0,ii), we have fti{t) € (0, it) and 
Pi(t) > for all i. Moving z m to the origin, we can define the (positive) cone created by the edges of 
the n-gon z m _iz m and z m z m+1 , as {a(z m _i - z m ) + b(z m+ i - z m ) : a,b > 0}. The n - 3 vertices 
which are not involved in creating the cone must lie in the interior of this cone for all t £ [0, ti). This 
is shown in Fig. By continuity of the Zj's, at time ti the vertices must lie either in the interior or on 
the boundary of this cone. But we have f3 m (h) = 0, implying that z m -\z m and z m z m +i are collinear 
and the cone is a line. Hence all the vertices are collinear, a contradiction by Lemma |6] ■ 
A straightforward consequence of the previous theorem is the following. 

Corollary 10: Consider an n-gon which is convex at t = 0. If the vertices evolve according to 
then for any t > 0, the n-gon will be strictly convex. 

Proof: Consider a vertex m for which /3 m (0) = it, and thus H m (0) = 0. We can choose this vertex 
such that H, m+ i(0) > since if Hi(0) = 0, Vi then the n-gon is not initially convex. From the proof 
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Z m -2> 



cone{z m _iz m , z m z m+ i} 
iZm+2 




Fig. 12. The evolution of a convex n-gon. The solid lines show the trajectories of each vertex. 



of Theorem [5] we have H m (t) = —2H m (t) + G m (t). But H m (0) = and we have shown in © that 
G m (0) > 0. Therefore, H m (0) > 0. By continuity of H m there exists a to > such that H m (t) > for 
f G [0,t ]- Thus, H m (t) > 0, for all t G (0,t ] and by Theorem^ fl" m (t) > for all t > t . ■ 
Fig. ^] shows the evolution of an initially convex n-gon. 



V. Optimal control law for perimeter shortening 

In [17] it is stated that a curve evolving according to Q is shrinking as fast as it can using only local 
information. To see why and in what sense this is true, reparametrize the curve in terms of its Euclidean 
arc-length s, defined via the differential arc-length element ds := \\dx/dp\\dp. With this we can write 
the length of a curve as 

Ut) r \ 9x 



Lit) 



ds 



Op 



dp. 



(9) 



In what follows we will differentiate this expression and determine the direction for the curve evolution 
which maximizes the rate of decrease of L(t). In order to take the time derivative of this expression, 
first consider taking the time derivative of ||cbc/<9p|| 2 : 



d_ 

di 



<9x 



dp 



d I cbc <9x 
dt\dp' dp 



I <9x d cbc 
\dp' dp dt 
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where (• , •) is the inner product (for u, v € W 1 , (u,v) = u T v). We also have that 

d_ 

at 

Therefore, combining these expressions and using the notation ||<9x/c?p|| = ||x||, we have 



<9x 


2 

= 2 


<9x 


-( 


<9x 


dp 




dp 


dt V 


dp 



d_ 

dt 

Using this expression in Q we obtain 

dL 
~dt 

Now, 



d 



x 



dp 



1 / <9x d d~> 



x|| \dp' dp dt I ' 



1 1 / cbc 3 <9x \ 



o 



x|| \dp' dp dt 



1 ds <9x cbc 



1 dx 



|x|| dp ||x|| 9p 9s ds ' 



since (is = ||x|| tip. This gives us 



Integrating by parts we obtain 

dL 
~dt 



dL 
~dt 

<9x <9x 
ds' , ~dt 







1 / <9x 9 9x\ 
\ds J dp dt / 

1 ^ <9 /<9x\ T <9x 
dp \ds J dt 



dp. 



o 

The first term on the right-hand side is zero since the curve is smooth and x(0,t) = x(l, i), and the 
second term can be rewritten to obtain 

T 



dL 



L 'd 2 x 



ds 2 



<9x 

Tt ds - 



Finally, since dx/ds = T and dT /ds = fcN, we have 

dL rL 
~dt ~ 



(10) 



Therefore, the direction of dx/dt in which L(t) is decreasing most rapidly is dx/dt = fcN, which is the 
Euclidean curve shortening rule ©. Note that this flow is optimal only in the sense that, given a fixed 
magnitude of the velocity of the curve at each point, this velocity always points in the direction which 
maximizes the rate of decrease of L(t). 

We now give an analogous result for the discrete polygon case. Given an n-gon we can write its 
perimeter as 



(ii) 



i=i 



In order to take the time derivative of P(t) consider taking the derivative of \zi + \ — Zi\ 2 
Zi,Zi + i — (for u, v € C n , (u, v) = u*v, where * denotes complex conjugate transpose). This yields 



d 

dt 



\ z i+l — Zi 



d_ 

dt 



(z-i+l — Zi, Zi+l — Zi) 



But also, 



dt 



\ z i+l — Zi 



2^.{{z i+ i - z h z i+1 - z^} . 
d 



Letting = U{ for % = 1, . . . , n and rearranging we have 
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Fig. 13. Evolving in the optimal direction. 



Therefore 



i=l 



p(t)=y^\( . Zi+1 Zi . , Ui+1 

\Zi+l — Zi\ 



Since all indices are evaluated modulo n this can be rewritten as 



n 



f-j' {\\Zi-i-Zi\ \Zi+x-Zi\ I) 

To maximize the rate of decrease of P(t), Ui should point in the direction of {z,i-\ — Zi)/\zi-\ — Zi\ + 
(zi+l — Zi)/\zi + i — Zi\. This direction bisects the internal angle 0i of the n-gon. In general, neither the 
linear scheme © nor the shortening by Menger-Melnikov curvature points in this direction. However, 
this direction does not ensure that the polygon becomes circular (nor elliptical); in simulation, adjacent 
vertices may capture each other and the polygon may collapse to a line.An example is shown in Fig. El 
Notice that these undesirable features do not appear in the linear scheme. 

Using dT2l and © we can determine P(t). For P(t) to be defined we require that adjacent vertices 
be distinct. This is ensured, for example, if the vertices start in a star formation about their centroid. The 
following result is analogous to the result in [13] that under Q the length of the curve monotonically 
decreases. 

Theorem 11: Consider an n-gon whose distinct vertices evolve according to ®. If adjacent vertices 
remain distinct, the perimeter P(t) of the n-gon monotonically decreases to zero. 
Proof: Substituting © into (fT^t and expanding we obtain 



Zi 



H+l 



Zi+l 



\ z i+l ~ z i\ 

Each term in this summation has the form — \u\ — \v\ + (u/\u\,v) + (v/\v\,u)}. From the Cauchy- 
Schwarz inequality we have ?R{(u/\u\,v}} < \v\,3l{{v/\v\,u)} < \u\, and thus \u\ — \v\ + (u/\u\, v) + 
(v/\v\,u)} < 0. Therefore, P(t) < 0. Equality is achieved if and only if u/\u\ = v/\v\ for each term in 
the summation; that is, if and only if 

Zi ~_ Zi ~\ = Zi+1 ~ Zt V*. (13) 

\Zi Zi—i\ \Zi-\-l Zi | 
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(a) The evolution of a simple polygon. The dashed 
lines show the trajectories of the vertices. 



(b) A plot of the area as a function of time. Notice 
that the area is initially increasing. 



Fig. 14. An embedded polygon for which the area initially increases. 



However, assume by way of contradiction that (THTl is satisfied. Rotate the coordinate system such that 
Z\ and Z2 lie on the real axis and Z2 — z\ > 0. Setting % = 2 in dl3l > we have z% — Z2 > 0, setting i = 3 
we have z^ — z% > 0, and so on. Hence m+i — Z{ > 0, Vi = 1, . . . , n — 1, which implies that z n > z\. But 
setting i = n in (fT3t we have z\ — z n > 0, a contradiction. Therefore (fT3t cannot be satisfied, P(t) < 0, 
and since the vertices converge to their stationary centroid, P(t) monotonically decreases to zero. □ 



There are two ways in which the linear scheme does not mimic Euclidean curve shortening. First of 
all, if an embedded curve is evolved via Euclidean curve shortening, its area is monotonically decreasing. 
However, for the linear scheme, in general, the area of a simple polygon is not monotonically decreasing. 
This is shown in Fig. ^] An interesting observation is that if a convex polygon evolves according to ©, 
its area is monotonically decreasing. To see this, consider a convex polygon at time t = with vertices 
Zi, i = 1, . . . , n, evolving according to ©. For each i, ij(0) is either zero, or points into the interior of 
the polygon, with ij(0) ^ for some i. Therefore, the area is initially decreasing. By Corollary ITOl the 
polygon is strictly convex for all t > 0, and thus points into the interior of the polygon for all i 
and for all t > 0. Therefore, the area decreases for all time. 

The second way in which the linear scheme does not mimic Euclidean curve shortening is in its effect 
on simple n-gons. If an embedded curve evolves according to the Euclidean curve shortening flow, it 
remains embedded. In contrast, a simple n-gon can become self-intersecting under the linear scheme, 
as is shown in Fig. ^] However, this is to be expected since the vertices in Fig. |15(a)| are not equally 
spaced around the polygon. The regions of the polygon with smaller spacing between adjacent vertices 
will move more slowly than the regions where the spacing is large. This is why, in Fig. |15(b)| the outer 
edge of the boomerang has intersected the inner edge. 



In summary, under the simple distributed linear control law (O, the robots rendezvous and also become 
more organized, in the sense that the polygon becomes elliptical. Furthermore, star formations remain 
so, convex polygons remain so, and the perimeter of the polygon decreases monotonically. These results 
are intended as a possible starting point for more useful behavior. As an example scenario, consider a 
number of mobile robots initially placed at random, and which should self-organize into a regular polygon 



VI. Limitations of the linear scheme 



VII. Conclusion 
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(a) A simple polygon. The vertices are marked by (b) The initial polygon evolves to the self-intersecting 

*'s. polygon shown by the thick solid line. 

Fig. 15. A simple polygon becomes self-intersecting. 



(circle) for the purpose of forming a large-aperture antenna. Distributed control laws for this goal would 
have to be nonlinear. Research on this front is on-going. 

Another topic for future research is to look at polygon shortening flows for wheeled robots which are 
subject to nonholonomic motion constraints. 

Finally, drawing upon the results on curve shortening flows, there has been a similar development of 
curve expanding flows — If a smooth, closed, and embedded curve is deformed along its outer normal 
vector field at a rate proportional to the inverse of its curvature, it expands to infinity, and the limiting 
shape is circular [24]. Thus, a scheme for deployment of a fleet of mobile robots could be achieved by 
creating the analogous polygon expanding flow. 
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